Dispersive Stabilization 

Guy Metivier? Jeffrey Rauch 1, 



Abstract 

111 posed linear and nonlinear initial value problems may be stabi- 
lized, that it converted to to well posed initial value problems, by the 
addition of purely nonscalar linear dispersive terms. This is a stabil- 
ity analog of the Turing instability. This idea applies to systems of 
quasilinear Schrodinger equations from nonlinear optics. 



1 Introduction 

In nonlinear optics, one commonly encounters coupled systems of scalar 
Schrodinger equations 

N 

(1.1) d tUj + iXjA xUj = b i,k(u, d x )u k , j e {1, . . . , N}, (t, x) G R 1+d , 

k=l 

where the \j are real and the bj ik are first order partial differential opera- 
tors with coefficients depending smoothly on u (see [2] and the references 
therein). The nonlinear terms usually depend on u and u, 

N 

(1.2) d t Uj + i\jA x Uj = c j>k (u, d x )u k + d jtk (u, d x )u k , 

k=l 

where the Cj tk and dj k are first order in d x . Introducing u and u as unknowns 
reduces to the form (j 1 . 1 j) for a doubled real system. 

For the local in time existence of smooth solutions, the easy case is when 
the first order part, B(u,d x )u on the right hand side is symmetric. In this 
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symmetric case there are easy L 2 estimates, followed by H s estimates ob- 
tained by commutations, which imply the local well-posedness of the Cauchy 
problem for (jl.ip in Sobolev spaces H s (R d ) for s > 1 + |. 

In many applications, B(u,d x ) is not symmetric and even more dt — 
B(u, d x ) is not hyperbolic and the Cauchy problem for dtu — B(u, d x )u = 
can be as ill posed as the Cauchy problem for the Laplacian. However, the 
Cauchy problem for (jl.ip may be well posed even if it is ill posed for the 
first order part. This is so even though the dispersive terms i\jA are not at 
all dissipative. We call this phenomenon dispersive stabilization. 

Example 1.1. With x G R the Cauchy problem for the system, 

o t u + i— -7j + d x v = 0, d t v - 1—7: - O x u = 0, 

is well posed in H s even though the first order part defines a badly ill posed 
initial value problem. This is proved by Fourier transformation in x. The 
amplification matrix is 

'if 

-ie 

For large ^ the matrix has purely imaginary eigenvalues close to ±i^ 2 and is 
uniformly diagonalisable showing that the amplification matrix is uniformly 
bounded for £ G R and t belonging to compact sets. The bound grows 
exponentially in time. The growth comes from |£| < R. 

The fact that the addition of a term diag (id 2 , — id x ) whose evolution 
is neutrally stable can stabilize a stongly ill posed Cauchy problem is not 
intuitively clear. There are many related results of this sort. The simplest is 
the following assertion about linear constant coefficient ordinary differential 
equaitons in the plane. 

Example 1.2. If A and B are 2x2 real matrices, knowing the stability 
origin as equilibrium of 

X' = AX, and, X' = B X, 

one can draw no conclusion about the stability of the equilibrium X' = 
(A + B)X. The best know is the Turing instability |15j for which A and B 
have eigenvalues with strictly negative real part so the input dynamics are 
exponentially stable and the sum dynamics can be unstable. Each of the 
stable dynamics is dissipative for certain scalar products. When the scalar 
products are different the Turing instability is possible. One but not both 
of the matrices A, B can be symmetric. 
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A related example is the two dimensional wave equation. 
Example 1.3. For the system version of the 2 — d wave equation, 

eac/i o/ i/ie sp/ii dynamics 

u t + (J U* = 0, ut + ^ % = 

defines a bounded semigroup on L°°(M 2 ). The first (resp. second) conserves 

\\uiWl 00 , and ||u2||l°°, ^resp. \\m + «2||i°°5 and \\u\ — ^([l 00 ) ■ 
The sum defines a dynamics so that the map 

u(0,x,y) i ^ u(t,x,y) 
is unbounded on L°°(M?) for allt^O. 

This analysis in this paper resembles example 11.11 We do not use the 
local smoothing properties of Schrodinger equations. Instead, the Fourier 
transform method is extended using the paradifferential calculus. The idea 
is to conjugate %A — B by a change of variable I + V with V of order —1 to 
a normal form 

(1.3) (ld + V){iA- B)(ld + Vy 1 =%A-B 

up to zero-th order terms, with B = i[V, A] — B symmetric. The conjugation 
(jl.3p means that the principal symbols satisfy 



(1.4) erg = a B +i[a A ,oy\. 

Equivalently, the energy estimates are obtained using the pseudodifferential 
symmetrizers 

(1.5) S = Id + V* + V 

If the Xj are pairwise distinct, one can reduce B to its diagonal part to prove 
the following result. 
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Theorem 1.4. If the Xj are real and pairwise distinct and if the diago- 
nal terms bjj(u,d x ) have real coeficients, then locally in time, the Cauchy 
problem for (jl.ip is well posed in the Sobolev spaces H s (R d ) for s > 1 + £. 

An analogous result for the systems (|1.2|) is the following. 

Theorem 1.5. Suppose that 

- the Xj are real and pairwise distinct 

- the diagonal terms bjj(u,d x ) have real coeficients, 

- Cj^Ujdx) = Ckj{u,d x ) for all pairs (J, k) such that Xj + = 0. 
Then locally in time, the Cauchy problem for (|1.2|) is well posed in the 
Sobolev spaces H s (R d ) with s > 1 + f . 

In the next section we give a more general statement which allows for 
more general nondiagonal second order terms. In particular the AjA^ can 
be replaced by different second order elliptic operators Aj(d x ). The idea 
of using pseudodifferential symmetrizers is related to the proof in [2] where 
the symmetry is obtained after differentiation of the equations and clever 
linear recombination. This amounts to using differential symmetrizers. Our 
analysis is a systematic exploration of the idea. Because of the quasilinear 
character of the equations, we use the paradifferential calculus in place of 
the classical pseudodifferential version. The latter would have sufficed to 
treat semilinear analogues. The paradifferential methods can also be used 
to treat the stongly nonlinear case F(u, d x u) since such a term is reduced 
to a quasilinear term using the paralinearization, see Section [21 

For the systems case the dispersive terms rotating at different speeds 
regularize an explosive first order term. For the scalar case, that is N = 1, 
such a stabilisation is not possible. The Cauchy problem for dt — iA x + id Xl 
is ill posed. However, if Im6(a;) satisfies suitable decay assumptions at in- 
finity, then the Cauchy problem for dt — i/S. x + b(x) ■ V x is well posed (see 
|12j). Intuitively, the waves propagate to the regions where b is small and 
are no longer amplified. The proofs use the dispersive and local smoothing 
properties of Schrbdinger equations. This idea has been extensively studied. 
Some of the foundational papers are [E], [6], [1], [7], [8], and, references 
therein. It would be natural to combine such ideas with those of disper- 
sive stabilization with the goal of extending the local existence to the case 
where the antisymmetric part of B has suitable decay at infinity rather than 
requiring that it vanish. We do not pursue this line of inquiry. 



4 



2 Statement of the result 

Consider the general equations, 

(2.1) d t u + iA(d x )u + B(t, x, u, d x )u = 0, 

with A second order and B first order, 

d 

(2.2) A{d x ) = A i,kd Xj d Xk , 

j,k=i 

d 

(2.3) B(t,x,u,d x ) = ^2Bj(t,x,u)d Xj . 

i=i 

The matrices Bj(t, x, u) are assumed to be C°° functions of (i, x, Rett, Im«), 
so that for each a and bouded K C C , 

d^uMuB G L°°([0,T\ xR d xK). 

With the example (11.11) in mind, we assume that A is smoothly block- 
diagonalizable: 

Assumption 2.1. For all £ G M n \{0}, A(^) = Yl^j,kCjCk is self-adjoint. 
Moreover, there are smooth real eigenvalues A p (£) and smooth self-adjoint 
eigenprojectors n p (£) such that 

^(0 = EM0n P (0- 

V 

This assumption is satisfied if A is self-adjoint with eigenvalues of con- 
stant multiplicity. The assumption allows for regular crossing of eigenvalues. 
The conditions on B involve, 

ImB := j\B-B*). 

Assumption 2.2. For all p and q, 

a p (0 = \(0 => n p (0im£(t,x,u,on g (0 = o. 

In addition, there are smooth matrix valued functions V pq (t,x,u,£) so that 

(2.4) n p (Olm£(t,x,u,£)ri g (0 = (\ p (0-X q (0)V p Jt,x,u,0- 
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Remark 2.3. The condition (|2,4p is automatic where A p (£) ^ A 9 (£), it 
defines V^, j(? . Assumption (|2.2p contains two types of information. 

• For any £, if A is an eigenvalue of A(£) and II(£) the spectral pro- 
jector, then Tl(£)B(t,x,u,£)IL(£) is self adjoint. If the eigenvalue remains of 
constant multiplicity for £ near £, nothing more needs to be added for this 
polarization. In particular, if all the distinct eigenvalues A p (£) of A(£) have 
constant multiplicity, the Assumption 12.21 reduces to the condition that the 
matrices TLp(£)B(t, x, u, £)n p (£) are self-adjoint. 

• If the eigenvalue A splits into several eigenvalues A p (£) for £ near £, 
the condition (|2.4p means that not only n p (£)Im-B(i, x, u, £)n g (£) vanishes 
at £ and on the variety {A p = A g }, but also that A p (£) — A 9 (£) is a divisor. 
In particular, if II(£) denotes the spectral projector on the invariant space 
associated to the eigenvalues close to A, this condition is locally satisfied 
with V Pt q = whenever II(£) B(t, x, u, £) II(£) is self-adjoint. This is so since 



Remark 2.4. There is no assumption on the spectrum of B(t,x,u,£). 
In particular, dt + B may be nonhyperbolic and thus strongly unstable in 
Hadamard's sense. The dispersive term A has a stabilizing effect, provided 
that the condition in Assumption 12.21 is satisfied. For this reason models of 
this type appear often in the descriptions of instabilities, for example that of 
Raman. The dispersive stabilisaton regularizes to a well posed causal model 
albeit with the possibility of growth for moderate wave numbers as in the 
example. 

We show that under the Assumptions 12.11 and 12.21 the Cauchy problem 
for (|2.ip is well posed in H s for s > | + 1, locally in time. 

Theorem 2.5. If Assumptions \2.1\ and \2.^ hold, s > and, h G H s (R d ), 

there is T > and a unique solution u G C°([0, T];H S (R d )) of (|2TPp with 
u\ t =o = h. 

Example 2.6 (From [2j). A is block diagonal A = diag{A p Id p } with real 
A p (£) homogeneous of degree two and A p (£) ^ A 9 (£) for p ^ q and £ ^ 0. The 
second assumption is trivially satisfied if the diagonal blocks B p ' p vanish. 

For applications, it is interesting to make explicit the assumptions when 
the first order part depends on u, 

(2.5) dtu + iA(d x )u + B(t, x, u, d x )u + C(t, x, u, d x )u = 




so, n p im5n 9 = n p nimBnn 5 = o. 
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Introducing v = u as a variable and setting U = t (u,v), the equation reads: 

(2.6) d t U + iA(d x )U + B(t,x,u,d x )U = 

with 

(A{d x ) \ /B C> 

(2-7) A=\ , B= 

\ -4(3,)/ \c 

In this context, the Assumption 12.21 becomes the following. 

Assumption 2.7. For allp andq, TL p (£)Im.B(t,x,u,£)Tlq(£) vanishes when 
^p(0 = and U p (£)(C(t,x,u,£) — t C(t, x, u, £))ll g (£) vanishes when 

A p (£) + A 9 (£) = 0. In addition, there are smooth matrices V p-q (t,x,u,^) and 
W P) q(t,x,u,£) such that 

(2.8) U p (ImB)U q = (A p - X q )V m , 

(2.9) U p (C- t C)U q = (X p + X q )W p , q . 



Theorem 2.8. Under Assumptions \2.1\ and \2.1\ for s > | + 1 and h € 

i? s (R d ), taere is T > and a unique solution u G C°([0,T];H s (R d )) of 
(|2.5p wrata data wU=q = 



We briefly discuss the case of equations with fully nonlinear right hand 
side, 

(2.10) d t u + iA{d x )u + F(t, x, u, d x u) = 0, 

where F(t, x,u,v\, . . . , v^) is a smooth function of (i, x, Rett, Imu) and of 
(Rei?i, . . . ,Im.Vd). Our analysis relies on a paralinearization of the first 
order term, so that the analogues of B and C are 

(2.11) B(t, x, u, v,C) = ^2 ^d^fo x > u ' v ^ 

i 3 
dF 

(2.12) C{t,x,u,v,0 = V^— (t,x,u,v) 

^— ' OVj 



with 



did id did id 



dvj 2 dKe Vj 2 dim Vj dvj 2 <9Re Vj 2 <9Im Vj 

as usual. The new condition is that (|2.8p (|2.9p are satisfied with smooth 

matrices x, «, u) and W Pi<i {t, x, u, v). In this case, the Cauchy problem 

d 
2 



is well posed in H s for s > i + 2 
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3 Basic L 2 estimate 

We solve (|2.ip by Picard iteration. Consider first the linear problem, 

(3.1) d t u + iA(d x )u + B(t,x,a,d x )u = f, u\ t=0 = h, 
where 

(3.2) aeC w ([0,T};H s (R d )), d t a G C° w ([0, T\; H s - 2 (R d )) 

with s > | + 1 and C° ([0, T]; F CT ) denotes the space of functions which are 
continuous from [0, T] to equipped with the weak topology. 

Theorem 3.1. There are functions Cq and C\ so that the solution of (|3.ip 
satisfies 

(3.3) |K*)|| L2 ^^(^^^(ll^ll^ + ^ll/COll^') 
with 

( 3 - 4 ) Kq := ||a||xoo([ ,r|xM< J )> 

( 3 - 5 ) ^1 := 1 1 a 1 1 £°° ( [0,T] ;H S (K^)) + ll^ a llL oo ([0,T];_H" s - 2 (IR d ))- 

Lemma 3.2 (Conjugation). For |£| large, there is a smooth invertible matrix 
V—i(t, x, u, £), homogeneous of degree —1 in sitc/i f/iai 

(3.6) #(M,«,£) - [K-ifta^u, 0,^(0] 
is self adjoint and homogeneous of degree 1 in £. 
Proof. Set 

y -! :=E r^r np(Imi?)np 

so that 5 - [V-!,A] = Yl ^pBU p is self adjoint. □ 

Proof of Theorem \3.1[ Use the paradifferential calculus and the notations 
of Section 5. 

a) For simplicity denote by Bj(t,x) the matrix Bj(t,x,a(t,x)) and by 
B = B(t,x,a(t,x),£) the symbol J2£,jBj. Because s > 1 + |, (|3.2|) implies 
that £j G C°([0,T];H s ), d t Bj G C°([0, T]; H s - 2 ) and 

(3.7) l|- B jlll,«'([0,T];H«(K' i )) + ll^t- B illL°°([0,T];H s - 2 (R d )) ^ Cl(-^l)- 
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In particular, as a symbol, B belongs to the class introduced in Defi- 
nition [57TTJ Using the paralinearization Proposition 15.81 we see that f± := 
B(t, x, d x )u — Tisu satisfies 

(3-8) \\fi(t)\\i? <Ci(^i)K*)||l», 

and u satisfies the paralinearized equation: 

(3.9) d t u + iA(d x )u + T iB u = / + u\ t=0 = h, 

b) Similarly, use the simplified notation V(t, x, £) = V_i(t, x, a(t, x), £)C(£) 
where £ £ C 00 (M <i ) vanishes near the origin and is equal to 1 for |£| > 1. 
Note that V € T^ 1 and that for all a there are functions Co iQ and C\ a such 
that for all t € [0, T] and £ G M d . 

(3.10) ||^(t,',0Lcc <Co,«(tf )(i + |£|) l ' i - a| 

(3.11) ||0£^(t, -,0L.- 2 < Ci,a(^i)(i + l^l) 1 ^ 1 • 

Use a symmetrizer, 

(3.12) E := Id + fTV-^TVr+TU-A^)- 1 . 

By Proposition 15.21 and Remark 15.71 there is a constant Cq(Kq) which de- 
pends only on Kq such that 

\\T v u(t)\\ H1 <C (K )\\u(t)\\ L2 . 

Therefore, 

(En,n) i2 > ||u|| £2 - 2C (K ) \\u\\ L2 \\u\Ijj_j + l\\u\\ 2 H -i- 
Choose 7 = j(Kq) so that 

(3.13) {Eu,u) L2 > -IMli 2 - 
Then, with another constant Cq(Kq), 

(3.14) H S ^)|| L 2 < Co(^o)|h(i)|| L2 - 

c) Compute 

(3.15) j t (E{t)u{t),u(t)) L2 = 2Re(Zd t u,u) L2 + ([$,E]u,u) £a . 



9 



By Lemma I5T91 and Proposition [57T2J [d t , S] = [c?f , Ty] + [9t, Ty]* is bounded 
from I? to L 2 and, 

(3.16) [[d t ,E]u(t),u(t)) L2 < d^Ollu^H^ . 

Next, observe that TyA^) = A(D X )T V + [Ty, A(a x )] = -T VA ^y There- 
fore, the equation and the symbolic calculus of Proposition 15.51 imply that 

m t u = -A{3 x )u - iA(d x )T v A{d x ) + i(T v )*A(D x ) - iT § u + £/ + f 2 

where B is the symbol B(t,x,£) — [V(t, x, £), A(£)] £ T\ and f 2 satisfies an 
estimate similar to (|3.8p . By Lemma 13.21 B is self adjoint for |£| > 2, and 
hence Proposition 15.61 implies that 

Re (iT B u(t),u(t)) L2 < CiTOK*)!^ • 

Since A(d x ) is self adjoint, we conclude that 

(3.17) j f [V(t)u(t),u(t)) L , < 2\\Xf(t)\\ L2 \\u(t)\\ L2 + d (KJWu^Wl, . 
Equations f)3. 13|) and (|3.14p imply estimate (|3.3|) . □ 



4 Sobolev estimates and nonlinear existence 



H s estimates for the linearized equation (|3.ip are obtained by differenti- 
ating the equation. The commutators [d x , B(t, x, a, d x ]u are estimated by 
standard nonlinear estimates as in the analysis of first order hyperbolic 
equations. Because s > i + 1, for \a\ < s, one has, 



(4.1) 



B(t,x,a,d x )]u(t) ^ <Ci(ifi)||«(t)| 



This implies the following estimates. 



Proposition 4.1. There are functions Cq and C\ such that smooth solutions 
of (|3.ip satisfy 

(4.2) < C (K )e tc ^ (\\u(0)\\ Hs + J*\\f(t')\\ Hs dt> 



with Kq and K\ defined at (|3.4p and (|3.5p . 



As in the hyperbolic theory, this estimates implies the following strong 
continuity result. 
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Proposition 4.2. Suppose that a satisfies (|3.2p , / G L 1 ([0, T], and 
h G F s . J/u G C°([0,T];iF) is a solution of ([33]) . tknuG C°([0, T], ff s ). 

Proof. With J e = (1 — eAj;) -1 , one checks that J £ u satisfies 

(4.3) d t J £ u + iA(d x )J e u + B(t,x,a,d x )J £ u = f £ , J e u\ t=0 = J e h, 

with / £ — > / in L 1 ([0, T], ff s ). Applying the estimates to J e u, one ob- 
tains that J £ u is a Cauchy family in C°([0, T], fP), implying that u G 
C°([0,T],H s ). □ 

Turn to the proof of the main result. More details can be found in |10j . 

Proof of Theorem \2.h\ (i) To solve f)3.1|) for a satisfying (|3.2[) use the molli- 
fied equations 

(4.4) d t u £ + iA{d x )J e u £ + B(t,x,a,d x )J £ u £ = f, u\ t=0 = h, 

where J e = (1 — eA. x )~ 1 . For fixed e, this is a linear o.d.e in since 
A(D X )J £ and i?J e are bounded. One checks that the proof of the esti- 
mates fj4.2|) for the solutions of (|3.ip immediately extends to the solutions 
of (|4.4p . because {J £ } is a bounded family of pseudodifferential operators 
of degree 0, are the new commutators they generate are remainders in the 
symbolic calculus developed in section [3l Therefore, the u £ are uniformly 
bounded in C°([0, T]; H s ). The equation shows that they are bounded in 

Extracting a subsequence and passing to the weak limit yields a solution 
u G C°([0,T],# s ). By Proposition H21 u G C°([0, T], ^ s ). 

(ii) Solve the nonlinear equation using the iteration scheme, 

(4.5) d t u n+ i + iA(d x )u n+ i + B(t,x,u n d x )u n+1 = 0, u n+1 \ t=0 = h. 

Using the estimate (I4.2p . one proves that there is T > such that the 
sequence {u n } is bounded in C°([0, T],H S ) and in C 1 ([0, T],H S ~ 2 ). Know- 
ing this bound in high norm, one checks that the sequence u n converges 
in a low norm C°([0,T];L 2 ). Passing to the limit gives a solution of (|2.ip 
tt G C°([0,T],F s ), which also belongs to C 1 ([0, T], fF~ 2 ). Using Proposi- 
tion 021 one obtains that u G C°([0,T],ff s ). □ 

5 Handbook of paradifferential calculus 

The symmetrizers are paradifferential operators in the variables x, depend- 
ing on the parameter t. This section reviews the paradifferential calculus 
extended to the case of time dependent symbols. 
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5.1 The spatial calculus 

Consider operators on M. d . The variables are denoted x and the frequency 
variables £. 

Definition 5.1 (Symbols). Let /i£l. 

i) Tq denotes the space of locally L°° functions a(x,^) on R d x R d which 
are C°° with respect to £ and such that for all a € N d there is a constant 
C a such that 

(5.1) V(x,0, \d?a(x,0\ < C«(l + |£|)"-H. 

iij r^ 1 denotes the space of symbols a £ such that for all j, d Xj a £ Tq. 

The paradifferential calculus in M. d , was introduced by J. M. Bony pQ (see 
also [11], [5], [H], [9]). The reference [10] gives a detailed account of the 
time dependent results needed here. The calculus associates operators T a to 
symbols a € Tq. They act in the scale of Sobolev spaces H S (W 1 ). Moreover, 
there is a symbolic calculus at order one for symbols in r^. Recall here the 
definition, as we will need it later on. 

Consider a C°° function ip(r),£) on M n x M n such that 

1) there are e± and 82 such that < E\ < 82 < 1 and 

f52l U(V,0 =1 for \ V \ <ei(l + |£|) 

\^(V,0 =0 for \rj\ >e 2 (l + |e|) • 

2) for all (a,/3) 6 N n x N n , there is C7 Qj/3 such that 

(5.3) v(»/,e,7) : i^f^^e^)! < c^ii + . 

For instance one can consider with iV > 3: 

+00 

k=0 

where x G C£°(R d ) satisfies < x < 1 and 

(5.5) *(£)=! for |^| < 1.1, X (0=0 for |^| > 1.9, 
and for k E Z, 

(5.6) x fc (e) = x(2- fc e) 
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and 

(5.7) ip = xo and for k > 1 tp k = Xk - Xk-i- 

A function ip satisfying (|5.2p (|5,3p is an admissible cut-off. Consider next 
■ , £) the inverse Fourier transform of • , £). For a € Tq define 

(5.8) := J G^(x-y,0a(y,0dy 
or equivalently on the Fourier side in x, 

(5-9) af(v,0 = Hv,Oa(r),0. 

The symbol a G Tq and belongs to Hormander's class S^r ^he P ar adiffer- 
ential operator T<f is defined by 

(5.10) lfu(x) := -±-J e**o*{x,S)u{(i)dZ. 

We collect here the main results. 

Proposition 5.2 (Action). Suppose that tp is an admissible cut-off. 

i) When a(£) is a symbol independent of x, the operator Ta is equal to 
the Fourier multiplier a(D). 

ii) For all a £ Tq and s £l, T<f is a bounded operator from H s (M. d ) to 

Proposition 5.3. If ip\ and ip% are two admissible cut-off, then for all a € 
rft and s£R, Tf 1 -Tf 2 is a bounded operator from H s (R d ) to H s ~^ +1 (R d ). 

Remark 5.4. This proposition implies that the choice of ijj is essentially 
irrelevant in our analysis, as in pQ. To simplify notation, make a definite 
choice of tp, for instance i/j = tpN with A = 3 as in (j5.4[) and use the notation 
T a for Tt 

Proposition 5.5 (Symbolic calculus). Consider a € and b € Tj 4 . Then 
ah € r^" 1- ^ and /or a// s € R, T a o — T a j, is bounded from H s (R d ) to 

7/ 6 is independent of x, then T a o Tj = T a b . 
These results extend to matrix valued symbols and operators. 
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Proposition 5.6 (Adjoints). Consider a matrix valued symbol a € r^. 
Denote by (T a )* the adjoint operator of T a in L 2 (M. d ) and by a*(x,£) the 
adjoint of the matrix a(x,£). Then (T a )* — T a * is bounded from H s (M. d ) to 



Remark 5.7. The norm of the operators acting in the indicated Sobolev 
spaces are uniformly bounded when the symbols a and b belong to bounded 
subsets of the symbol classes. 

Bounded functions of x are particular examples of symbols in the class 
rf), independent of the frequency variables (. In this case, T a is called a 
paraproduct in pQ. 

Proposition 5.8 (Paralinearization). There is a constant C such that for 
all a e W^°° and all u G l?( 



\\ad Xj u - T a d Xj u\\ L2 < C\\a\\wi,oo ||w||i,2 . 
5.2 The time dependent case 

In the sequel we consider functions of (£, x) € [0, T] x M ra , considered as 
functions of t with values in various spaces of functions of x. In particular, 
denote by T a the operator acting on u so that for each fixed t, (T a u)(t) = 
T a (t)u{t). 

(5.11) T a u(t,x) := J e**<r a (t,x,£)u(0dZ. 
with 

(5.12) a a (t,x,0 := J G(x-y,€)a(t,y,Z)dy 

This definition shows that formally 

(5.13) [d t ,T a ] = T dta . 

This yields easy estimates when dta £ L°° . However, if we want to keep the 
lower bound s > 1 + | in Theorem 12.51 this condition need not be satisfied, 
since in the equation (|2.ip . dt has the weight of two spatial derivatives. This 
is why we introduce a slight extension. 

Using the Littlewood-Paley decomposition 

+oo 

(5.14) u = AfcM, with Aku:=(fkU, 

k=0 
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as in (|5.7p . the Besov space Boo '°° is denned as the space of tempered 
distributions u such that 

(5.15) ||u|| R -i,cx; = sup 2 _fc ||Afcn|L 00 < +00. 

II II D^Q II I I Li 

This space occurs in our analysis because of the the following embedding. 

Lemma 5.9. Functions u £ H s belong to i?^, 1 ' 00 when s > i — 1. 

In the spirit of Definition 15-H introduce the following notation. 

Definition 5.10. For fi £ M, letT l ^ 1 denote the space of distributions a(x , £) 
on M. d x W 1 which are C°° with respect to £ with values in B^' 00 and such 
that for all a £ N d there is a constant C a such that 

(5-16) V£, \\d^a(;0\\ B -i,oo < C a (1 + leir -1011 . 

Definition 5.11 (Time dependent symbols). Let fj, £ R and T > 0. 

i) Tq denotes the space of locally continuous functions a(t, x, £) on [0, T] x 
M rf xR rf which are C°° with respect to £ and such that the family {a(t, • , ■ );t £ 
[0, T]} is bounded in Tq. 

denotes the space of symbols a £ Tq suc/i i/iarf 

- i/ie family {a(t, • , • ); i £ [0, T]} is bounded in 

- the family {dta(t, ■ , ■ );t £ [0, T]} is bounded in . 

For a £ T^, the operator T a is defined by (|5.1ip and the Propositions 15.21 
15.51 15.81 apply for fixed t, yielding estimates that are uniform in t (see Re- 
mark [52]). The commutation with d% is treated as follows. 

Proposition 5.12. For a £ , the commutator [dt, T a ] maps C°([0, T]; H s ) 
to C°([0,T];H s -» +1 ) and there is a constant C such that for all t € [0,T] 

(5.17) \\[B t ,T a ]u{t)\\ H .^. 1 <C\\u\\ H '. 

Moreover, the constant C depends only on the supremum for t £ [0, T] of a 
finite number of semi-norms 

(5.18) sup(l + |£|)l a Mto(-,£)|L-i jOC . 

^ IIS II Dqo 

Proof. One has, 

+00 

(5.19) d t a a (t,-,0 = ^ S k - N (D x )(dta(t,-,S))<p k (0. 

k=0 
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where Sj is the Fourier multiplier with symbol Xj- Since dta(t, •,£) belongs 

to Boo 1 ' 00 , 

\S k _ N (D x ){d t a(t, < 2 fe (i + ieir. 

On the support of ip k , that is |£| ~ 2 k . With similar estimates for the 
derivatives, this shows that dta(t,-,-) is a bounded family of symbols in 
S^ 1 . By construction, the spectral property that dta a {t,ViO 1S supported 
in \r)\ < e(l + |£|) for some e > is satisfied and therefore the operator 
(dta)(t,x, D x ) is bounded from H s to H s ~^~ 1 for all s. □ 
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